Systematics of threshold incident energy for deep sub-barrier fusion hindrance 
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We systematically evaluate the potential energy at the touching configuration for heavy-ion reactions using 
various potential models. We point out that the energy at the touching point, especially that estimated with the 
Krappe-Nix-Sierk (KNS) potential, strongly correlates with the threshold incident energy for steep fall-off of 
fusion cross sections observed recently for several systems at extremely low energies. This clearly indicates that 
the steep fall-off phenomenon can be attributed to the dynamics after the target and projectile touch with each 
other, e.g., the tunneling process and the nuclear saturation property in the overlap region. 

PACS numbers: 25.60.Pj, 24.10.Eq, 25.70Jj,25.70.-z 



Recently, for medium-heavy mass systems, it has become 
possible to measure fusion cross sections down to extremely 
low incident energies. In those measurements, unexpected 
steep fall-off of fusion cross sections, as compared to a stan- 
dard theoretical calculation, have been observed at deep sub- 
barrier energies U H, H, 0, 13, Si • Although the steep fall-off 
phenomenon, referred to as the fusion hindrance, may be ac- 
counted for if one uses an anomalously large diffuseness pa- 
rameter in the Woods-Saxon potential |2[, the physical origin 
of the phenomenon has yet to be clarified [8]. 

One important aspect of fusion reactions at deep subbarrier 
energies is that the inner turning point of the potential may be 
located far inside the touching point of the colliding nuclei. 
We show this schematically in Fig. 1. At energies close to the 
Coulomb barrier, the inner turning point is still far outside of 
the touching point ||9[] (see the line (i) in Fig.QJ. At these en- 
ergies, one usually assumes that a compound nucleus is auto- 
matically formed once the projectile penetrates the Coulomb 
barrier, due to the strong nuclear attractive force in the clas- 
sically allowed region. In contrast, at energies below the po- 
tential energy at the touching point, Vrouch, the inner turning 
point appears more inside of the touching point (see the line 
(ii) in Fig.[T|). That is, the projectile nucleus is still in the clas- 
sically forbidden region when the two colliding nuclei touch 
with each other. After the touching, an elongated composite 
system is formed, which evolves in the classically forbidden 
region towards a compound nucleus by overlapping between 
the projectile-like and the target-like fragments. Since this 
involves the penetration of the residual Coulomb barrier, nat- 
urally the fusion cross sections are hindered by the tunneling 
factor. 

In this paper, we evaluate the potential energy at the touch- 
ing configuration for several systems, and investigate whether 
the dynamics after the touching point is responsible for the 
steep fall-off phenomenon. In this respect, it is interesting 
to notice that the authors of Refs. [1,4] have argued that the 
steep fall-off phenomenon systematically takes place below 
a certain threshold incident energy, E s . We will show below 
that there is a strong correlation between the touching energy 
^Touch an d the threshold energy E s , indicating that the den- 
sity overlap in the classically forbidden region indeed plays 
an important role. We mention that one would have to settle 



i— 
01 

c 
LU 



c 

s 

o 

Q_ 



Overlap Region 
Sudden \ 



Compound ?' 
State S 

» ' Adiabatic 



Two-Body Region 




Touching 
Point 



_ 



Center-of-Mass Distance r 

FIG. 1: (Color online) Schematic picture for heavy-ion sub-barrier 
fusion reactions. The filled circle denotes the energy at the touching 
point, Vrouch- 



a model in the overlap region, such as the adiabatic or sudden 
models, or some combination of these two, in order to clarify 
the whole dynamics of deep subbarrier fusion reactions. How- 
ever, our analysis is independent of these modellings, since 
both the adiabatic and the sudden approaches provide a simi- 
lar potential energy to each other as long as the touching point 
is concerned. 

In order to estimate the potential energy at the touching 
point, ?Wch = Rp + Rr, we employ the Krappe-Nix-Sierk 
(KNS) ini, the Bass lH, the proximity (U and the Akyiiz- 
Winther (AW) lfl3ll models. Assuming the spherical shape for 
both the projectile and target nuclei, the KNS potential energy 
at the touching point reads 



V W - —D 

V KNS — U 



rtouch fiRr/a) f(Rp/a) 
a g{Rr/a) g(Rp/a) 



(1) 



1 sinh(jc) 



where the functions / and g are defined as f(x) 
and g(x) — xcosh(;t:) — sinh(;t:), respectively. In this model, the 
nuclear radius is given by R = roA 1 / 3 , and the depth constant 
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TABLE I: Potential energy at the touching configuration calculated by various theoretical models for the systems discussed in Ref. [5|]. In 
the second column, ZjZp is the charge product of the system, while fl is the reduced mass. X e ff is the effective fissility parameter defined 



as X e {{ = ZpZp/[A T A p (A p +A T )]/12. 6, while E s is the experimental energy at which the astrophysical S-factor has the maximum [3] 
^KNS> ^Prox> ^Bass an d V"aw denote the results of the Krappe-Nix-Sierk, the proximity, the Bass and the Akyuz-Winther models, respectively. 
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0.705 


175 ± 1.8 
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0.331 
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27.7 


21.9 


35.4 
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47.3 ± 0.9 
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0.282 


27.6 ± 0.8 


26.1 
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16 Q+ 16q 
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7.1 ±0.8 
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where the effective surface energy constant c s is given by 
c s = a s (l - K S I 2 ) with I = (N — Z)/A. We take the param- 
eters to be a = 0.68 fm, a s = 21 .33 MeV and K s = 2.378 from 
FRLDM2002 lfl2ll . except for the radius parameter for which 
we slightly adjust to be rrj = 1.2 fm in order to fit the ex- 
perimental fusion cross sections for the 64 Ni+ 64 Ni reaction at 
energies above the Coulomb barrier. 

The proximity potential energy at the touching point is 
given by 



yW = _ 1 . 7818 ^_3.00 ) 



(3) 



where R — RtRp/(Rt +Rp)- In this model, the nuclear radius 
is given by 

/ 7 b 2 49 b 4 \ N 



with 



Roo = 1.240A 1/3 [1 + 1.646/A-0.191(A-2Z)/A)], (5) 



3 J(N-Z)/A-±c 1 ZA- 1 / 3 
' ~ 2 r ° 6+|7A-V3 



(6) 



The value of the parameters are taken to be b=l fm, ro = 1 . 14 
fm, 7=32.65 MeV, ci=0.757895 MeV and g=35.4 MeV. The 
surface energy coefficient 02 in Eq. ([3]) is given by 02 = 
18.36 — Q(tj + tp) /2rp. In order to fit the experimental data, 
we use the same prescription as in Refs. il5lll6ll and subtract 
3.00 MeV from the original proximity model (the last term in 
Eq. ®). 

The Bass potential energy at the touching point is given by 



■P] 



-1 



(7) 



where the parameters a and j3 are taken as a=0.0300 MeV -1 
fm and j3=0.0061 MeV -1 fm, respectively. In the Bass model, 
the nuclear radius is given bytf= 1.16A l / 3 -139A- 1 / 3 . The 
AW potential energy at the touching point, on the other hand, 
reads 



Vj§ = -SnyRa, 



(8) 



where the average surface tension 7 is given by 7 = 



0.95 



1-1. 



}( N T -Z T \i Np-Zp ' 



MeV fm . In this model, 



the nuclear radius is given by R = (l^OA^-O^) fm 
and the diffuseness parameter a is given by a = 0.855 • [1 + 



O^A^+A^ 1 '")]- 1 fm. 

In order to estimate the total potential energy at the touch- 
ing point, Vrouch. one nas to add the Coulomb potential to 
the nuclear potential energies given by Eqs. (fj) - ®. To 
this end, we use the Coulomb potential for two point charges, 
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FIG. 2: (Color online) Potential energy at the touching point cal- 
culated by various theoretical models. The open circles, triangles, 
diamonds and squares denote the results of the Krappe-Nix-Sierk, 
the Bass, the Proximity and the Akyiiz-Winther models, respectively. 
The solid line denotes the systematics proposed by Jiang et al. [5j]. 
The filled circles, squares, triangles, and the horizontal lines show the 
experimental energy taken from Ref . 1 5] at which the astrophysical 
S-factor has the maximum value. 



V ( c ) = e 2 ZjZp /rt uch, where the touching radius r t0U ch is spec- 
ified for each model for the nuclear potential, V^>. The re- 
sultant touching energy for the systems discussed in Ref. [5] 
is shown in Fig.[T]as a function of ZrZpji 1 ! 2 , where jX is the 
reduced mass of the colliding nuclei. The results of the KNS, 
the Bass, the proximity and the AW models are denoted by the 
open circle, the open triangle, the open diamond and the open 
square, respectively. All the results are summarized in Table 
1. These touching energies are compared with the energy E s , 
at which the experimental fusion cross section is maximum 
when it is plotted in terms of the astrophysical S-factor yj]. 
These "experimental" energies E s are shown in Fig. 2 by the 
filled circles, the filled squares, the filled triangles, and the 
horizontal lines, depending on the types of the system as de- 
fined in Ref. |5|]. Notice that the energy E s for the type III was 
estimated by extrapolation, and that for the type IV is only an 
upper limit. The systematics for the energy E s proposed by 
Jiang et al. is also shown by the solid line. 

Althogh the physical significance for the energy E s is not 
clear, because the S-factor representation for fusion cross sec- 
tions would be useful only at much lower energies than the 
lowest energies of the current measurements, at which the 
outer turning point is much larger than the inner turning point 



TABLE II: The potential energy at the touching configuration calcu- 
lated with the KNS model for the systems discussed in Refs. (2,0,0]. 
All of these systems are categolized as Type III. 
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1704. 
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\/.JOJ 


4Q 




r?i 
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28 + 62jvjj 


1722 


^70 


48 (\ 


44 1 


r?i 
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16q + 208p b 


2529 


0.413 


69.6 


70.5 


[3] 


16 0+ 144 Sm 


1882 


0.384 


57.7 


54.6 


[3] 


19p + 208p b 


3079 


0.431 


75.5 


78.8 


[3] 


4(, Ca + 90 Zr 


4210 


0.524 


93.2 


93.6 


[3] 


50 Ti + 208p b 


11454 


0.683 


181.2 


191.9 


[3] 


58 Ni + 60 Ni 


4258 


0.527 


92 ±2 


92.5 


[4] 


58 Ni + 64 Ni 


4325 


0.511 


89 ±2 


91.1 


[4] 


58 Ni + 74 Ge 


5109 


0.542 


98.5 ±2.0 


103.9 


[4] 


64 Ni + 74 Ge 


5249 


0.517 


97.5 ± 2.0 


101.5 


[4J 



(see e.g., Refs. C3, H for a discussion on the modified S- 
factor, that takes into account the effect of the inner turning 
point), it is remarkable that the result of the KNS model fol- 
lows closely to the energy E s , and thus the systematics shown 
by the solid line (an exceptional case of 64 Ni+ 100 Mo will be 
discussed in the next paragraph). The good correspondence 
between Vkns and E s may be due to the fact that the KNS 
model partly takes into account the saturation of nuclear mat- 
ter when two nuclei come inside the Coulomb barrier fll9ll (in 
fact, the KNS model has been shown to be consistent with 
the ener gy d ensity formalism with the Skyrme SkM* interac- 
tion 1201 Bill ). The result of the AW potential is similar to that 
of the KNS model, although the deviation from E s is slightly 
larger. For the Bass and the proximity models, although the 
dependence of the touching energy Vxouch on the parameter 
ZjZp^l 1 is similar to that of the KNS and the AW models, 
there is a large discrepancy between the touching energy Vt OU ch 
and the threshold energy E s . 

For the asymmetric 64 Ni+ 100 Mo reaction, the experimen- 
tal threshold energy E s deviates largely from the systematics 
curve. The calculations with the KNS and AW models are 
consistent with the systematics curve but not to the value of 
E s . In order to check how the touching energy VVouch compares 
with the threshold energy E s for other asymmetric systems, we 
also examine the 16 O+ 208 Pb reaction. For this system, we find 
that the KNS model leads to the touching energy Vxouch that 
is consistent with the experimental threshold energy E s JH 
(see Table II). Therefore, it is unlikely that the large differ- 
ence between Vxouch and E s for the 64 Ni+ 100 Mo system can 
be attributed to the model assumption of the KNS potential. 
Notice that for this system, there may exist some peculiar nu- 
clear structure effect, because the coupled-channels calcula- 
tion reported in Ref. [22] does not seem to account well for 
the experimental fusion cross sections even above the thresh- 
old energy E s . A further investigation is necessary for this 
system concerning the threshold energy. 

In order to see more clearly the correlation between E s 
and Vkns, the lower panel of Fig. [3] shows these energies as 
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FIG. 3: (Color online) Same as Fig. 2, but as a function of the effec- 
tive fissility parameter defined in Ref. 1 1 1]. The upper panel shows 
the difference between E s and Vkns ■ 



that the difference between E s and Vkns is indeed small ex- 
cept for large X e g, clearly indicating that there is a strong cor- 
relation between these two values. The large discrepancy for 
systems with large X e ff may be due to the ambiguity of the ex- 
perimental data, because the measurements were only for the 
fusion-evaporation cross sections and the fusion-fission cross 
sections were estimated using the statistical model JUS 51]. 

In summary, we have shown that the potential energy at 
the touching point strongly correlates with the threshold in- 
cident energies for the steep fall-off of the fusion cross sec- 
tions. The systematics of the threshold energy can be rather 
naturally explained by the present approach in terms of the 
touching energy. This strongly suggests that the overlap pro- 
cess after the touching is responsible for the steep fall-off of 
the fusion cross section. For such overlap process, the sudden 
and adiabatic approaches have been often employed 11231 124I1 . 
In the former, the frozen density approximation while over- 
lapping with the colliding nuclei is applied, and in the latter, 
the dynamical change in the density of the colliding nuclei is 
taken into account. These two approaches are in the opposite 
limit to each other, and there is not yet a definite consensus 
regarding which limit better describes the realistic situation at 
deep subbarrier energies. In this respect, the threshold energy 
discussed in this paper will provide a useful constraint to mod- 
elling of the overlap process as the touching configuration is a 
doorway of such process. 



a function of the effective fissility parameter X e g defined as 



X, 



cff 



ZTZp/lA^Ap^iA]/ 3 



Aj/ 3 )]/12.6 |H1. The figure 
includes also a few more systems than shown in Fig. 2, which 
are taken from Refs. |Q> H Si (see Table II for the additional 
data). With this representation, all the data points distribute 
more uniformly than in Fig. [2] The difference between E s and 
Vkns is also shown in the upper panel of Fig. [3] One observes 
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